3.1.66 \(\int \frac {x^3 \sin (c+d x)}{(a+b x^2)^2} \, dx\) [66]

3.1.66.1 Optimal result

Integrand size = 19, antiderivative size = 431 \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx=\frac {\sqrt {-a} d \cos \left (c+\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{4 b^{5/2}}-\frac {\sqrt {-a} d \cos \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}+d x\right )}{4 b^{5/2}}+\frac {\operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}+d x\right ) \sin \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b^2}+\frac {\operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right ) \sin \left (c+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b^2}+\frac {\sin (c+d x)}{2 b^2}-\frac {x^2 \sin (c+d x)}{2 b \left (a+b x^2\right )}-\frac {\cos \left (c+\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b^2}+\frac {\sqrt {-a} d \sin \left (c+\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{4 b^{5/2}}+\frac {\cos \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}+d x\right )}{2 b^2}+\frac {\sqrt {-a} d \sin \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}+d x\right )}{4 b^{5/2}} \]

output

1/2*cos(c+d*(-a)^(1/2)/b^(1/2))*Si(d*x-d*(-a)^(1/2)/b^(1/2))/b^2+1/2*cos(c 
-d*(-a)^(1/2)/b^(1/2))*Si(d*x+d*(-a)^(1/2)/b^(1/2))/b^2+1/2*sin(d*x+c)/b^2 
-1/2*x^2*sin(d*x+c)/b/(b*x^2+a)+1/2*Ci(d*x+d*(-a)^(1/2)/b^(1/2))*sin(c-d*( 
-a)^(1/2)/b^(1/2))/b^2+1/2*Ci(-d*x+d*(-a)^(1/2)/b^(1/2))*sin(c+d*(-a)^(1/2 
)/b^(1/2))/b^2-1/4*d*Ci(d*x+d*(-a)^(1/2)/b^(1/2))*cos(c-d*(-a)^(1/2)/b^(1/ 
2))*(-a)^(1/2)/b^(5/2)+1/4*d*Ci(-d*x+d*(-a)^(1/2)/b^(1/2))*cos(c+d*(-a)^(1 
/2)/b^(1/2))*(-a)^(1/2)/b^(5/2)+1/4*d*Si(d*x+d*(-a)^(1/2)/b^(1/2))*sin(c-d 
*(-a)^(1/2)/b^(1/2))*(-a)^(1/2)/b^(5/2)-1/4*d*Si(d*x-d*(-a)^(1/2)/b^(1/2)) 
*sin(c+d*(-a)^(1/2)/b^(1/2))*(-a)^(1/2)/b^(5/2)
 

3.1.66.2 Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 1.82 (sec) , antiderivative size = 286, normalized size of antiderivative = 0.66 \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx=\frac {i e^{-i c-\frac {\sqrt {a} d}{\sqrt {b}}} \left (\left (2 \sqrt {b}+\sqrt {a} d\right ) e^{\frac {2 \sqrt {a} d}{\sqrt {b}}} \operatorname {ExpIntegralEi}\left (-\frac {\sqrt {a} d}{\sqrt {b}}-i d x\right )+\left (2 \sqrt {b}-\sqrt {a} d\right ) \operatorname {ExpIntegralEi}\left (\frac {\sqrt {a} d}{\sqrt {b}}-i d x\right )\right )-i e^{i c-\frac {\sqrt {a} d}{\sqrt {b}}} \left (\left (2 \sqrt {b}+\sqrt {a} d\right ) e^{\frac {2 \sqrt {a} d}{\sqrt {b}}} \operatorname {ExpIntegralEi}\left (-\frac {\sqrt {a} d}{\sqrt {b}}+i d x\right )+\left (2 \sqrt {b}-\sqrt {a} d\right ) \operatorname {ExpIntegralEi}\left (\frac {\sqrt {a} d}{\sqrt {b}}+i d x\right )\right )+\frac {4 a \sqrt {b} \cos (d x) \sin (c)}{a+b x^2}+\frac {4 a \sqrt {b} \cos (c) \sin (d x)}{a+b x^2}}{8 b^{5/2}} \]

input

Integrate[(x^3*Sin[c + d*x])/(a + b*x^2)^2,x]
 

output

(I*E^((-I)*c - (Sqrt[a]*d)/Sqrt[b])*((2*Sqrt[b] + Sqrt[a]*d)*E^((2*Sqrt[a] 
*d)/Sqrt[b])*ExpIntegralEi[-((Sqrt[a]*d)/Sqrt[b]) - I*d*x] + (2*Sqrt[b] - 
Sqrt[a]*d)*ExpIntegralEi[(Sqrt[a]*d)/Sqrt[b] - I*d*x]) - I*E^(I*c - (Sqrt[ 
a]*d)/Sqrt[b])*((2*Sqrt[b] + Sqrt[a]*d)*E^((2*Sqrt[a]*d)/Sqrt[b])*ExpInteg 
ralEi[-((Sqrt[a]*d)/Sqrt[b]) + I*d*x] + (2*Sqrt[b] - Sqrt[a]*d)*ExpIntegra 
lEi[(Sqrt[a]*d)/Sqrt[b] + I*d*x]) + (4*a*Sqrt[b]*Cos[d*x]*Sin[c])/(a + b*x 
^2) + (4*a*Sqrt[b]*Cos[c]*Sin[d*x])/(a + b*x^2))/(8*b^(5/2))
 

3.1.66.3 Rubi [A] (verified)

Time = 1.00 (sec) , antiderivative size = 441, normalized size of antiderivative = 1.02, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.263, Rules used = {3824, 3826, 2009, 3827, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx\)

\(\Big \downarrow \) 3824

\(\displaystyle \frac {\int \frac {x \sin (c+d x)}{b x^2+a}dx}{b}+\frac {d \int \frac {x^2 \cos (c+d x)}{b x^2+a}dx}{2 b}-\frac {x^2 \sin (c+d x)}{2 b \left (a+b x^2\right )}\)

\(\Big \downarrow \) 3826

\(\displaystyle \frac {d \int \frac {x^2 \cos (c+d x)}{b x^2+a}dx}{2 b}+\frac {\int \left (\frac {\sin (c+d x)}{2 \sqrt {b} \left (\sqrt {b} x+\sqrt {-a}\right )}-\frac {\sin (c+d x)}{2 \sqrt {b} \left (\sqrt {-a}-\sqrt {b} x\right )}\right )dx}{b}-\frac {x^2 \sin (c+d x)}{2 b \left (a+b x^2\right )}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {d \int \frac {x^2 \cos (c+d x)}{b x^2+a}dx}{2 b}+\frac {\frac {\sin \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b}+\frac {\sin \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b}-\frac {\cos \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b}+\frac {\cos \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b}}{b}-\frac {x^2 \sin (c+d x)}{2 b \left (a+b x^2\right )}\)

\(\Big \downarrow \) 3827

\(\displaystyle \frac {d \int \left (\frac {\cos (c+d x)}{b}-\frac {a \cos (c+d x)}{b \left (b x^2+a\right )}\right )dx}{2 b}+\frac {\frac {\sin \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b}+\frac {\sin \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b}-\frac {\cos \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b}+\frac {\cos \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b}}{b}-\frac {x^2 \sin (c+d x)}{2 b \left (a+b x^2\right )}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {d \left (\frac {\sqrt {-a} \cos \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b^{3/2}}-\frac {\sqrt {-a} \cos \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b^{3/2}}+\frac {\sqrt {-a} \sin \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b^{3/2}}+\frac {\sqrt {-a} \sin \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b^{3/2}}+\frac {\sin (c+d x)}{b d}\right )}{2 b}+\frac {\frac {\sin \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b}+\frac {\sin \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b}-\frac {\cos \left (\frac {\sqrt {-a} d}{\sqrt {b}}+c\right ) \text {Si}\left (\frac {\sqrt {-a} d}{\sqrt {b}}-d x\right )}{2 b}+\frac {\cos \left (c-\frac {\sqrt {-a} d}{\sqrt {b}}\right ) \text {Si}\left (x d+\frac {\sqrt {-a} d}{\sqrt {b}}\right )}{2 b}}{b}-\frac {x^2 \sin (c+d x)}{2 b \left (a+b x^2\right )}\)

input

Int[(x^3*Sin[c + d*x])/(a + b*x^2)^2,x]
 

output

-1/2*(x^2*Sin[c + d*x])/(b*(a + b*x^2)) + ((CosIntegral[(Sqrt[-a]*d)/Sqrt[ 
b] + d*x]*Sin[c - (Sqrt[-a]*d)/Sqrt[b]])/(2*b) + (CosIntegral[(Sqrt[-a]*d) 
/Sqrt[b] - d*x]*Sin[c + (Sqrt[-a]*d)/Sqrt[b]])/(2*b) - (Cos[c + (Sqrt[-a]* 
d)/Sqrt[b]]*SinIntegral[(Sqrt[-a]*d)/Sqrt[b] - d*x])/(2*b) + (Cos[c - (Sqr 
t[-a]*d)/Sqrt[b]]*SinIntegral[(Sqrt[-a]*d)/Sqrt[b] + d*x])/(2*b))/b + (d*( 
(Sqrt[-a]*Cos[c + (Sqrt[-a]*d)/Sqrt[b]]*CosIntegral[(Sqrt[-a]*d)/Sqrt[b] - 
 d*x])/(2*b^(3/2)) - (Sqrt[-a]*Cos[c - (Sqrt[-a]*d)/Sqrt[b]]*CosIntegral[( 
Sqrt[-a]*d)/Sqrt[b] + d*x])/(2*b^(3/2)) + Sin[c + d*x]/(b*d) + (Sqrt[-a]*S 
in[c + (Sqrt[-a]*d)/Sqrt[b]]*SinIntegral[(Sqrt[-a]*d)/Sqrt[b] - d*x])/(2*b 
^(3/2)) + (Sqrt[-a]*Sin[c - (Sqrt[-a]*d)/Sqrt[b]]*SinIntegral[(Sqrt[-a]*d) 
/Sqrt[b] + d*x])/(2*b^(3/2))))/(2*b)
 

3.1.66.3.1 Defintions of rubi rules used

rule 2009

Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 3824

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*Sin[(c_.) + (d_.)*(x_)], x_Sym 
bol] :> Simp[x^(m - n + 1)*(a + b*x^n)^(p + 1)*(Sin[c + d*x]/(b*n*(p + 1))) 
, x] + (-Simp[(m - n + 1)/(b*n*(p + 1))   Int[x^(m - n)*(a + b*x^n)^(p + 1) 
*Sin[c + d*x], x], x] - Simp[d/(b*n*(p + 1))   Int[x^(m - n + 1)*(a + b*x^n 
)^(p + 1)*Cos[c + d*x], x], x]) /; FreeQ[{a, b, c, d, m}, x] && ILtQ[p, -1] 
 && IGtQ[n, 0] && (GtQ[m - n + 1, 0] || GtQ[n, 2]) && RationalQ[m]
 

rule 3826

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*Sin[(c_.) + (d_.)*(x_)], x_Sym 
bol] :> Int[ExpandIntegrand[Sin[c + d*x], x^m*(a + b*x^n)^p, x], x] /; Free 
Q[{a, b, c, d, m}, x] && ILtQ[p, 0] && IGtQ[n, 0] && (EqQ[n, 2] || EqQ[p, - 
1]) && IntegerQ[m]
 

rule 3827

Int[Cos[(c_.) + (d_.)*(x_)]*(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Sym 
bol] :> Int[ExpandIntegrand[Cos[c + d*x], x^m*(a + b*x^n)^p, x], x] /; Free 
Q[{a, b, c, d, m}, x] && ILtQ[p, 0] && IGtQ[n, 0] && (EqQ[n, 2] || EqQ[p, - 
1]) && IntegerQ[m]
 

3.1.66.4 Maple [C] (verified)

Result contains complex when optimal does not.

Time = 0.63 (sec) , antiderivative size = 852, normalized size of antiderivative = 1.98

method	result	size
risch	\(\frac {i \sqrt {a b}\, {\mathrm e}^{\frac {i c b +d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (\frac {i c b +d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right ) d}{8 b^{3}}-\frac {i \sqrt {a b}\, {\mathrm e}^{\frac {i c b -d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (\frac {i c b -d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right ) d}{8 b^{3}}+\frac {i {\mathrm e}^{\frac {i c b +d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (\frac {i c b +d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right )}{4 b^{2}}+\frac {i {\mathrm e}^{\frac {i c b -d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (\frac {i c b -d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right )}{4 b^{2}}+\frac {i {\mathrm e}^{-\frac {i c b +d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (-\frac {i c b +d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right ) \sqrt {a b}\, d}{8 b^{3}}-\frac {i {\mathrm e}^{-\frac {i c b -d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (-\frac {i c b -d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right ) \sqrt {a b}\, d}{8 b^{3}}-\frac {i {\mathrm e}^{-\frac {i c b +d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (-\frac {i c b +d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right )}{4 b^{2}}-\frac {i {\mathrm e}^{-\frac {i c b -d \sqrt {a b}}{b}} \operatorname {Ei}_{1}\left (-\frac {i c b -d \sqrt {a b}-b \left (i d x +i c \right )}{b}\right )}{4 b^{2}}+\frac {\left (\frac {\left (3 i \left (i d x +i c \right ) a b c \,d^{2}-i \left (i d x +i c \right ) b^{2} c^{3}+a^{2} d^{4}-b^{2} c^{4}\right ) d^{2}}{2 \left (2 i \left (i d x +i c \right ) b c -b \left (i d x +i c \right )^{2}+a \,d^{2}+c^{2} b \right ) b^{2} a}+\frac {c^{3} d^{3} x}{2 \left (-2 i \left (i d x +i c \right ) b c +b \left (i d x +i c \right )^{2}-a \,d^{2}-c^{2} b \right ) a}-\frac {3 c^{2} d^{2} \left (i \left (i d x +i c \right ) b c +a \,d^{2}+c^{2} b \right )}{2 a b \left (2 i \left (i d x +i c \right ) b c -b \left (i d x +i c \right )^{2}+a \,d^{2}+c^{2} b \right )}-\frac {3 i c \,d^{2} \left (-i a c \,d^{2}-i b \,c^{3}-\left (i d x +i c \right ) a \,d^{2}+\left (i d x +i c \right ) b \,c^{2}\right )}{2 a b \left (-2 i \left (i d x +i c \right ) b c +b \left (i d x +i c \right )^{2}-a \,d^{2}-c^{2} b \right )}\right ) \sin \left (d x +c \right )}{d^{4}}\)	\(852\)
derivativedivides	\(\text {Expression too large to display}\)	\(2531\)
default	\(\text {Expression too large to display}\)	\(2531\)

input

int(x^3*sin(d*x+c)/(b*x^2+a)^2,x,method=_RETURNVERBOSE)
 

output

1/8*I/b^3*(a*b)^(1/2)*exp((I*c*b+d*(a*b)^(1/2))/b)*Ei(1,(I*c*b+d*(a*b)^(1/ 
2)-b*(I*d*x+I*c))/b)*d-1/8*I/b^3*(a*b)^(1/2)*exp((I*c*b-d*(a*b)^(1/2))/b)* 
Ei(1,(I*c*b-d*(a*b)^(1/2)-b*(I*d*x+I*c))/b)*d+1/4*I/b^2*exp((I*c*b+d*(a*b) 
^(1/2))/b)*Ei(1,(I*c*b+d*(a*b)^(1/2)-b*(I*d*x+I*c))/b)+1/4*I/b^2*exp((I*c* 
b-d*(a*b)^(1/2))/b)*Ei(1,(I*c*b-d*(a*b)^(1/2)-b*(I*d*x+I*c))/b)+1/8*I/b^3* 
exp(-(I*c*b+d*(a*b)^(1/2))/b)*Ei(1,-(I*c*b+d*(a*b)^(1/2)-b*(I*d*x+I*c))/b) 
*(a*b)^(1/2)*d-1/8*I/b^3*exp(-(I*c*b-d*(a*b)^(1/2))/b)*Ei(1,-(I*c*b-d*(a*b 
)^(1/2)-b*(I*d*x+I*c))/b)*(a*b)^(1/2)*d-1/4*I/b^2*exp(-(I*c*b+d*(a*b)^(1/2 
))/b)*Ei(1,-(I*c*b+d*(a*b)^(1/2)-b*(I*d*x+I*c))/b)-1/4*I/b^2*exp(-(I*c*b-d 
*(a*b)^(1/2))/b)*Ei(1,-(I*c*b-d*(a*b)^(1/2)-b*(I*d*x+I*c))/b)+1/d^4*(1/2/( 
2*I*(I*d*x+I*c)*b*c-b*(I*d*x+I*c)^2+a*d^2+c^2*b)/b^2/a*(3*I*(I*d*x+I*c)*a* 
b*c*d^2-I*(I*d*x+I*c)*b^2*c^3+a^2*d^4-b^2*c^4)*d^2+1/2*c^3*d^3*x/(-2*I*(I* 
d*x+I*c)*b*c+b*(I*d*x+I*c)^2-a*d^2-c^2*b)/a-3/2*c^2*d^2*(I*(I*d*x+I*c)*b*c 
+a*d^2+c^2*b)/a/b/(2*I*(I*d*x+I*c)*b*c-b*(I*d*x+I*c)^2+a*d^2+c^2*b)-3/2*I* 
c*d^2*(-I*a*c*d^2-I*b*c^3-(I*d*x+I*c)*a*d^2+(I*d*x+I*c)*b*c^2)/a/b/(-2*I*( 
I*d*x+I*c)*b*c+b*(I*d*x+I*c)^2-a*d^2-c^2*b))*sin(d*x+c)
 

3.1.66.5 Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.32 (sec) , antiderivative size = 291, normalized size of antiderivative = 0.68 \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx=-\frac {{\left (2 i \, b x^{2} - {\left (-i \, b x^{2} - i \, a\right )} \sqrt {\frac {a d^{2}}{b}} + 2 i \, a\right )} {\rm Ei}\left (i \, d x - \sqrt {\frac {a d^{2}}{b}}\right ) e^{\left (i \, c + \sqrt {\frac {a d^{2}}{b}}\right )} + {\left (2 i \, b x^{2} - {\left (i \, b x^{2} + i \, a\right )} \sqrt {\frac {a d^{2}}{b}} + 2 i \, a\right )} {\rm Ei}\left (i \, d x + \sqrt {\frac {a d^{2}}{b}}\right ) e^{\left (i \, c - \sqrt {\frac {a d^{2}}{b}}\right )} + {\left (-2 i \, b x^{2} - {\left (i \, b x^{2} + i \, a\right )} \sqrt {\frac {a d^{2}}{b}} - 2 i \, a\right )} {\rm Ei}\left (-i \, d x - \sqrt {\frac {a d^{2}}{b}}\right ) e^{\left (-i \, c + \sqrt {\frac {a d^{2}}{b}}\right )} + {\left (-2 i \, b x^{2} - {\left (-i \, b x^{2} - i \, a\right )} \sqrt {\frac {a d^{2}}{b}} - 2 i \, a\right )} {\rm Ei}\left (-i \, d x + \sqrt {\frac {a d^{2}}{b}}\right ) e^{\left (-i \, c - \sqrt {\frac {a d^{2}}{b}}\right )} - 4 \, a \sin \left (d x + c\right )}{8 \, {\left (b^{3} x^{2} + a b^{2}\right )}} \]

input

integrate(x^3*sin(d*x+c)/(b*x^2+a)^2,x, algorithm="fricas")
 

output

-1/8*((2*I*b*x^2 - (-I*b*x^2 - I*a)*sqrt(a*d^2/b) + 2*I*a)*Ei(I*d*x - sqrt 
(a*d^2/b))*e^(I*c + sqrt(a*d^2/b)) + (2*I*b*x^2 - (I*b*x^2 + I*a)*sqrt(a*d 
^2/b) + 2*I*a)*Ei(I*d*x + sqrt(a*d^2/b))*e^(I*c - sqrt(a*d^2/b)) + (-2*I*b 
*x^2 - (I*b*x^2 + I*a)*sqrt(a*d^2/b) - 2*I*a)*Ei(-I*d*x - sqrt(a*d^2/b))*e 
^(-I*c + sqrt(a*d^2/b)) + (-2*I*b*x^2 - (-I*b*x^2 - I*a)*sqrt(a*d^2/b) - 2 
*I*a)*Ei(-I*d*x + sqrt(a*d^2/b))*e^(-I*c - sqrt(a*d^2/b)) - 4*a*sin(d*x + 
c))/(b^3*x^2 + a*b^2)
 

3.1.66.6 Sympy [F]

\[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx=\int \frac {x^{3} \sin {\left (c + d x \right )}}{\left (a + b x^{2}\right )^{2}}\, dx \]

input

integrate(x**3*sin(d*x+c)/(b*x**2+a)**2,x)
 

output

Integral(x**3*sin(c + d*x)/(a + b*x**2)**2, x)
 

3.1.66.7 Maxima [F]

\[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx=\int { \frac {x^{3} \sin \left (d x + c\right )}{{\left (b x^{2} + a\right )}^{2}} \,d x } \]

input

integrate(x^3*sin(d*x+c)/(b*x^2+a)^2,x, algorithm="maxima")
 

output

-1/2*((cos(c)^2 + sin(c)^2)*d*x^2*sin(d*x + c) + ((d^2*x^3*cos(c) - d*x^2* 
sin(c) - 2*x*cos(c))*cos(d*x + c)^2 + (d^2*x^3*cos(c) - d*x^2*sin(c) - 2*x 
*cos(c))*sin(d*x + c)^2)*cos(d*x + 2*c) + ((cos(c)^2 + sin(c)^2)*d^2*x^3 - 
 2*(cos(c)^2 + sin(c)^2)*x)*cos(d*x + c) - 2*(((b^2*cos(c)^2 + b^2*sin(c)^ 
2)*d^3*x^4 + 2*(a*b*cos(c)^2 + a*b*sin(c)^2)*d^3*x^2 + (a^2*cos(c)^2 + a^2 
*sin(c)^2)*d^3)*cos(d*x + c)^2 + ((b^2*cos(c)^2 + b^2*sin(c)^2)*d^3*x^4 + 
2*(a*b*cos(c)^2 + a*b*sin(c)^2)*d^3*x^2 + (a^2*cos(c)^2 + a^2*sin(c)^2)*d^ 
3)*sin(d*x + c)^2)*integrate((2*a*d*x*sin(d*x + c) + ((2*a*d^2 + 3*b)*x^2 
- a)*cos(d*x + c))/(b^3*d^3*x^6 + 3*a*b^2*d^3*x^4 + 3*a^2*b*d^3*x^2 + a^3* 
d^3), x) - 2*(((b^2*cos(c)^2 + b^2*sin(c)^2)*d^3*x^4 + 2*(a*b*cos(c)^2 + a 
*b*sin(c)^2)*d^3*x^2 + (a^2*cos(c)^2 + a^2*sin(c)^2)*d^3)*cos(d*x + c)^2 + 
 ((b^2*cos(c)^2 + b^2*sin(c)^2)*d^3*x^4 + 2*(a*b*cos(c)^2 + a*b*sin(c)^2)* 
d^3*x^2 + (a^2*cos(c)^2 + a^2*sin(c)^2)*d^3)*sin(d*x + c)^2)*integrate((2* 
a*d*x*sin(d*x + c) + ((2*a*d^2 + 3*b)*x^2 - a)*cos(d*x + c))/((b^3*d^3*x^6 
 + 3*a*b^2*d^3*x^4 + 3*a^2*b*d^3*x^2 + a^3*d^3)*cos(d*x + c)^2 + (b^3*d^3* 
x^6 + 3*a*b^2*d^3*x^4 + 3*a^2*b*d^3*x^2 + a^3*d^3)*sin(d*x + c)^2), x) + ( 
(d^2*x^3*sin(c) + d*x^2*cos(c) - 2*x*sin(c))*cos(d*x + c)^2 + (d^2*x^3*sin 
(c) + d*x^2*cos(c) - 2*x*sin(c))*sin(d*x + c)^2)*sin(d*x + 2*c))/(((b^2*co 
s(c)^2 + b^2*sin(c)^2)*d^3*x^4 + 2*(a*b*cos(c)^2 + a*b*sin(c)^2)*d^3*x^2 + 
 (a^2*cos(c)^2 + a^2*sin(c)^2)*d^3)*cos(d*x + c)^2 + ((b^2*cos(c)^2 + b...
 

3.1.66.8 Giac [F]

\[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx=\int { \frac {x^{3} \sin \left (d x + c\right )}{{\left (b x^{2} + a\right )}^{2}} \,d x } \]

input

integrate(x^3*sin(d*x+c)/(b*x^2+a)^2,x, algorithm="giac")
 

output

integrate(x^3*sin(d*x + c)/(b*x^2 + a)^2, x)
 

3.1.66.9 Mupad [F(-1)]

Timed out. \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^2\right )^2} \, dx=\int \frac {x^3\,\sin \left (c+d\,x\right )}{{\left (b\,x^2+a\right )}^2} \,d x \]

input

int((x^3*sin(c + d*x))/(a + b*x^2)^2,x)
 

output

int((x^3*sin(c + d*x))/(a + b*x^2)^2, x)